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Abstract 


We perform an analysis of semi-inclusive production of charm (D-mesons) in neutral 
current (NC) and charged current (CC) deep inelastic scattering (DIS) at full 0(a]). Our 
calculation is based on the heavy quark scheme developed by Aivazis, Collins, Olness and 
Tung (ACOT) where we include an 0(a]) calculation of quark scattering contributions for 
general masses and couplings. We review the relevant massive formulae and subtraction 
terms and discuss their massless limits. We show how the charm fragmentation function 
can be measured in CC DIS and we investigate whether the charm production dynamics 
may be tested in NC DIS. We also discuss finite initial state quark mass effects in CC 
and NC DIS. 



1 Introduction 


In a recent paper Jl[] we have analyzed heavy quark initiated contributions to fully inclusive 
deep inelastic (DI) structure functions. Towards lower values of the Bjorken variable x, 
heavy (charm) quarks are produced in about 20 % of the neutral current (NC) li 
and charged current (CC) [Q deep inelastic events in lepton-nucleon collisions. Therefore 
in this kinematical range, heavy quark events contribute an important component to 
the fully inclusive DI structure functions of the nucleon. However, due to acceptance 
losses this component can usually not be measured directly by inclusively tagging on 
charm events and more differential observables have to be considered like E'd I, 1,0, 
Pt or p [|2|, |3j] distributions, where Ed, px and p are the energy, transverse momentum 
and pseudorapidity of the charmed hadron produced, i.e. mainly of a meson. In 
this article we consider Ed spectra represented by the usual scaling variable z defined 
below. Within DIS the charmed hadron energy spectrum is the distribution which is most 
sensitive to the charm fragmentation process and may give complementary information 
to one hadron inclusive e + e~ annihilation which is usually chosen to define fragmentation 
functions (FF’s) || |j. A well understanding of charm fragmentation is essential for any 
charm observable, e.g. the normalization of px and p distributions in photoproduction is 
substantially influenced by the hardness of the FF |10[ [lTJ . The z distribution of charm 
fragments is directly measured in CC neutrinoproduction 1&00 and may give insight 
into details of the charm production dynamics in NC electroproduction. It has e.g. been 
shown in |2j. [3|, that the energy spectrum of ZA*)-mesons produced at the eP-collider 
HERA may be able to discriminate between intrinsic and extrinsic production of charm 
quarks. 


Intrinsic heavy quark densities may arise due to a nonperturbative component of 
the nucleon wave function [|13|] or due to a perturbative resummation Jl4|, [15], 


| 17] of 


large quasi-collinear logs [ln(Q 2 /m 2 ); Q and m being the virtuality of the mediated gauge 
boson and the heavy quark mass, respectively] arising at any order in fixed order extrinsic 
production (or fixed order perturbation theory: FOPT). Here we will only consider the 
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latter possibility for inducing an intrinsic charm density c(x, Q 2 ) which is concentrated 
at small x and we will ignore nonperturbative components which are expected to be 
located at large x 0. Technically the resummation of large perturbative logs proceeds 
through calculating the boundary conditions for a transformation of the factorization 
scheme which is switched from rtf to rif + 1 active, massless flavors, canonically 

at Q 2 = m 2 . For fully inclusive DIS the Kinoshita-Lee-Nauenberg theorem confines all 
quasi-collinear logs to the initial state such that they may be absorbed into c(x, Q 2 ). For 
semi-inclusive DIS (SI DIS) also final state collinearities arise which are resummed in 
perturbative fragmentation functions D?(z, Q 2 ) (parton i decaying into charm quark c) 
along the lines of [|H| [HJ. The scale dependence of c(x, Q 2 ) and D^{z, Q 2 ) is governed by 


massless renormalization group (RG) evolution. 


Besides this zero mass variable flavor number scheme, where mass effects are only 
taken care of by the boundary conditions for c(x, Q 2 ) and Dflz, Q 2 ), variable flavor number 
schemes have been formulated [F|, [15|, [16|, [17]], which aim at resumming the quasi-collinear 
logs as outlined above while also keeping power suppressed terms of 0[{m 2 / Q 2 ) k ] in the 
perturbative coefficient functions. Our reference scheme for this type of schemes will be 
the one developed by Aivazis, Collins, Olness and Tung (ACOT) j2U|. |H| . In the ACOT 
scheme full dependence on the heavy quark mass is kept in graphs containing heavy 
quark lines. This gives rise to the above mentioned quasi-collinear logs and to the power 
suppressed terms. While the latter are regarded as mass corrections to the massless, 
dimensionally regularized, standard coefficient functions (e.g. in the MS scheme), the 
former are removed numerically by subtraction terms, which are obtained from the small 
mass limit of the massive coefficient functions. 


The outline of our paper will be the following: In section 2 we will shortly overview the 
relevant formulae for SI DIS for general masses and couplings including quark scattering 
(QS) and boson gluon fusion (GF) contributions up to 0(a]). We will thereby present our 
ACOT based calculation for the O,S'- 11 [] component of SI structure functions. In section 
1 Bracketed upper indices count powers of a s . 
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3 we will analyze the charm fragmentation function in CC and NC DIS. In section 4 we 
draw our conclusions and some uncomfortably long formulae are relegated to an appendix. 


2 Semi-Inclusive Heavy Quark Structure Functions 


This section presents the relevant formulae for one (heavy flavored) hadron inclusive DIS 
structure functions. The contributions from scattering events on massive quarks are given 
up to 0(a]) in section 2T and GF contributions are briefly recalled in section |2d^ . Section 
[2.3| presents all subtraction terms which render the structure functions infrared safe and 
includes a discussion of these terms. 


2.1 Scattering on Massive Quarks 


We consider DIS of the virtual Boson B* with momentum q on the quark Qi with mass 
mi and momentum pi producing the quark Q 2 with mass m 2 and momentum P 2 . The 
latter fragments into a heavy quark flavored hadron Hq 2 , e.g. a \Q 2 Q 1 ) meson, q L being any 
light quark. Phenomenologically most prominent are of course charm quarks fragmenting 
into Th*)-mesons which are the lightest heavy flavored hadrons. 


We will strictly take over the formulae and notations of our inclusive analysis in [4]] 
whenever possible and merely extend them for SI DIS considered here. In particular 
we take over the definition of the structure functions Ti given in terms of the usual 
experimental structure functions F) in Eq. (7) in |[J, i.e. 


Ti 


2A 

S' + E ++ — 2mim 2 S_ 

S + A 2x 2 



Fi 


( 1 ) 


with 


E ±± = Q 2 ± 7712 ± m\ 


( 2 ) 
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In Eq. ([!;) we use the shorthand A = A— Q 2 ] , where the usual triangle function 
is defined by 

A [a, 6, c] = \J a 2 + b 2 + c 2 — 2(afc + be + ca) . (3) 

The vector (V) and axial vector (A) couplings of the — Ajs)Qi quark current 

enter via the following combinations: 

S± = VV' ± AA' 

R± = (VA' ± V'A)/2 (4) 

where V, A = V', A' in the case of pure B scattering and V, A ^ V', A' in the case of B , B' 
interference (e.g. 7 , Z° interference in the standard model). 


Since we want to investigate the energy spectrum of charm fragments, we introduce 
the Lorentz-invariant z = Ph Q2 • Pn/q • Pn which reduces to the energy E Hq scaled to its 
maximal value v = q 0 in the target rest frame. Therefore in contrast to Eq. (A4) in [|T| we 
do not integrate the tensor over the full partonic phase space but keep it differential 
in the corresponding partonic variable z r = P2 ■ Pi/q ■ Pi 01 the mass corrected variable z 
which is defined below. In order to obtain hadronic observables we have to extend the 


ansatz of Eq. (17) in [20] such that it includes a nonperturbative hadronization function 
Dq 2 . In the limit of vanishing masses the massless parton model expressions have to be 
recovered. Our ansatz will be 

di d( 


= 


f C 


QiiZit *) d Q 2 (C, [p ]) |{ p + = ^p+. 2= ^| 


(5) 


where /i is the factorization scale, 5 = z' / z' LO with z' LO = E ++ /E + _ and v + = (v° +v 3 )/V2 
for a general vector v. W^ v is the usual hadronic tensor and is its partonic analogue. 
Eq. (|5|) defines the fragmentation function Dq 2 to be a multiplicative factor multiplying 
inclusive structure functions at LO/Born accuracy, i.e. 

vQ 2 ) = Q 2 ) d q&- B 2 ]) . 


(0 


where the fpCii (x,Q 2 ) are defined and given in |T|. The scale dependence of Dq 2 is 
bracketed here and in the following because it is optional; a more detailed discussion on 
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this point will be given at the end of section |2.3| . We do not construct onr ansatz in Eq. 
(U) for the convolution of the fragmentation function along light front components for the 
outgoing particles which would only be Lorentz-invariant for boosts along a specified axis. 
Since the final state of DIS is spread over the entire solid angle it has no preferred axis 
as defined for the initial state by the beam direction of collider experiments. Note that 
eq. (|5|) is in agreement with usual factorized expressions for massless initial state quanta, 
e.g. in 


2 lH , since there mi = 0 such that z' LO = 1. 


Up to 0(a>\) the hadronic structure functions for scattering on a heavy quark read 

ay(/r 2 ) f 1 f 1 dz 

! x Y L t 


U Q 2 ; = Qi(x, tf) A? 2 (u [/^D + 


2vr 


(7) 


x 


Q &n (X.c 1 ) 8, 


D 


Q 2 \ W 


with 20 


X = 


x 

2 Q 1 


+ A 


( 8 ) 


In Eq. (|T|) and throughout this paper we set the renormalization scale equal to the fac¬ 
torization scale. The kinematical boundaries of the phase space in Eq. ([7[) are introduced 
by the theta function cut Q q . In the massless limit Q q —> 1. The precise arguments of 0 q 
are set by the kinematical requirement 


z! ■ < z' < z' 

min max i 


(9) 


with 


Z max — 
min 


±A[s,ml, -Q 2 ](s - ml) + (Q 2 + m 2 + s)(s + m 2 ) 


2s(s — m\ + Q 2 ) 

where s = (pi + q) 2 . Note that Eq. (^) also poses an implicit constraint on via z' min < 


The Hf for nonzero masses are obtained in exactly the same way as outlined for fully 
inclusive structure functions in [|l| if the partonic phase space is not fully integrated over. 
They are given by 

H q Zi',z') = C F [(Si + K) 5(1-0 5(1-z) 
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( 10 ) 


1 — €' z LO Si + E 


_ j ;LO 

"" (1-0+ 8 

with the normalization factors 

S + £ ++ — c lm\m2S_ 


F 1 "r 1 


A 


N i = 


2A 


N o = 


and the mandelstam variables 


2S' + A 
A' ’ 


1V S = 


2R + 
“A 7 " ’ 


where A = 


£i + E ++ 

Si + E + _ 


si(0 = ( Pi + qf-ml 

= b|^[(A - E + _)f'+ A + E + _] 

= (Pi-Paf-ml 
— [si(^ 7 ) + S + -](^ - z o) 

is the would-be pole of the /-channel propagator and we use A' = 


A[s, mf, — Q 2 ]. Si and Vi are the soft real and virtual contribution to H respectively. 
They can be found in appendix C of [|lj whereas the fi(§i, ti) are listed in appendix A of 
this paper. In the massless limit the Hf reduce to the MS coefficient functions in B p2] 


up to some divergent subtraction terms which we will specify in section |273 . 


2.2 Gluon Fusion Contributions at 0(a]) 


The semi-inclusive coefficient functions for GF production of massive quarks have been 
obtained for general masses and couplings in |H|] and are only briefly reviewed here for 
completeness. They are given by 

F?J ( x , z, Q 2 , /i 2 ) = [ y g(x',fi 2 ) /i, 3 (|pC,<3 2 ) Dq 2 (^, [/r 2 ]) Q g 


F 2 gf (x,z,Q 2 ,jli 2 ) = 


dx’ f 1 d( , . , 


x 


^ x'g{x\p ) / 2 (—,C,Q ) Dq 2 {~, [h ]) 0 9 (11) 


with ax = [ 1 + (mi + m 2 ) 2 /Q 2 )x. The Q g cut guarantees that ( m i n < ( < (max where 
C min,max are given in Eq. (3) of J2lj]. Similarly to Eq. (||) (min < C max may also constrain 
the phase space available for the x' integration. 
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2.3 Subtraction Terms 


It requires three MS subtraction terms to render the double convolutions in Eqs. (0), (0) 
infrared safe: 




X 


V 

f 1 d? 

1 

+ 

to 


L e 

Li-e 

X 2 " 

Dq 2 ( 



In -0 — 1 — 2 ln(l — 

m\ 


srSUBg, 2 \ _ n ( r 2 n\ ^(Z^ 2 ) i ft 2 f d£ p( 0 )/£/\ X 5 

^ (x, 2, Q ? M ) -^Q2 [ft J ) 2 y|- ^ ^2 J a ^Q9 ) 9 ^ 7 ft 


( 12 ) 

(13) 


rr*-(x,z,M = 

x ^q 2 (~v 2 ) Qi{x,v 2 ) , 


L±^lAn4-l-21n(l-^) 

1 — Z' \ 777-2 


(14) 


where Pqg( £,') = 1/2 [£ /2 + (1 — £') 2 ]. Note that SUB g in Eq. (p~3|) differs slightly from Eq. 
(6) in [flip because we are allowing for a nonzero initial state parton mass m\ here which 


we did not in 21 


The subtraction terms define the running of the initial state quark density ( SUB q , 
SUBg) and the final state fragmentation function ( SUBd ) in the massless limit. They 
remove collincar logarithms and scheme defining finite terms from the convolutions in 
Eqs. (0), (|TT]) and they are constructed such that the massless MS results of P, |2|] are 
recovered in the limit 


lim 

mi,2 —>0 . 


r, 


QS^+V+GF 


0, z, Q 2 , fi 2 ) - T\ 


SUB q +SUBg+SUB D 


(x,z,Q 2 ,fi 2 ) 


^ 1)MS (x,z,Q 2 ) 


(15) 


where SUB q and SUB D regularize fi® s( + J whereas SUB g regularizes Tf F . Contrary to 
the fully inclusive SUB g term in Eq. (16) of Q there is no need to include an additional 
~ ln( / u 2 /??72) subtraction in Eq. ( |T3D because the </ _1 h-channel singularity of the massless 
limit of GFC) is located at £ = 0 and is outside the integration volume of Eq. (|TT|). 
If only initial state subtractions, i.e. SUB q and SUB g , are considered and final state 
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subtractions, i.e. SUBp, are not performed one reproduces, in the limit m\ —> 0 the 
results in [[L2| for producing a heavy quark from a light quark. Note that in this case the 
fragmentation function Dq 2 should be taken scale-independent, say of the Peterson form 
|23| . In the limit where also the final state quark mass m 2 approaches zero and the final 
state subtraction term SUB D is subtracted from the results in JT2[] the massless quark 
results in ^ are obtained and a running of Dq 2 is induced via a RG resummation 
of final state collinear logs as formulated for one hadron inclusive e + e _ annihilation in 

§§ 

Apart from removing the long distance physics from the coefficient functions the sub¬ 
traction terms set the boundary conditions for the intrinsic heavy quark density Q 1 [|l| 
and the perturbative part of the heavy quark fragmentation function Dq 2 |[K)| 


Qi{x,Qo) = 


2vr 


"■3 /'?«»(!•« 


(16) 


d Qi (z,qd = 2 #> C F r 1<u 


2vr 


1+z 


/ 2 


In 


Q. 


mt> 


- 1 — 2 ln(l - z!) 


D 


02 ' z' 


(17) 


where Qq, Qq are the transition scales at which the factorization scheme is switched from 
rif to rif + 1, rif+ 2 active flavors, respectively (assuming here for simplicity that m\ < m 2 ; 
a generalization to m\ > m 2 is straightforward). For general Qq also the gluon density and 
a s undergo a scheme transformation. Here we omit the corresponding formulae which can 
be found in [|R| [L7| since they are not needed in our analysis. Canonically Qq is set equal 
to the heavy quark mass rri\ which guarantees [[H| up to two loops a continuous evolution 
of a s and the light parton densities across Qq. All available heavy quark densities are 
generated using Q 0 = m\ and we will therefore follow this choice here although a variation 
of Qq might substantially influence the heavy quark results even far above the threshold 
m • Note that at three loops a continuous evolution across Qq can no longer be achieved, 
neither for the parton distributions 


nor for a s |25|. Analogously to Q 0 = m± we 
use Qq = m 2 throughout. In Eq. (0) we have made the distinction between the scale 
dependent FF Dq 2 (z, Qq) and the scale independent FF Dq 2 (z ) explicit. Following the 


























terminology in ||T9| , p~0[| the latter corresponds to the nonperturbative part of the former 
and describes the hadronization process at the end of the parton shower which is described 
perturbatively by the massless RG evolution. Alternatively, Dq 2 (z) corresponds to a scale- 
independent FF within FOPT where no collinear resummations are performed. These two 
points of view may induce a scheme dependence if Dq 2 ( z ) is fitted to data. In principle, 
the massless evolution equations generate nonzero FF’s also for light partons to decay 
into heavy flavored hadrons. These light—>heavy contributions are important at LEP 
energies [JTD[ but can be safely neglected at the scales considered hercQ and we will assume 
Di^Q 2 (z, /i 2 ) = 0 throughout. 


2.4 SI Structure Functions at 0{a\) 

In the next section we will consider three types of O(al) VFNS structure functions 
first two are constructed at full 0{ot\) 

Q5(o+i) +GF j-.SU B q +SUB g +[SUB D \ 

^ i F i ? 

where the inclusion or omission of the bracketed SUBp term corresponds to a running 
or scale-independent fragmentation function, respectively, as discussed in the previous 
section. It is somewhat unclear whether QS ^ contributions (and the corresponding 
subtractions) should be considered on the same perturbative level as GF^ l \ see [p]] for a 
more detailed discussion on that point. In the original formulation of the ACOT scheme 
|T4| QS ^ contributions are neglected at the level we are considering here and we therefore 
do also consider this option via the partial 0(a]) structure function 

jpQSW+GF _ pSUBg . 

For obtaining the numerical results of the next section the general formulae of this section 
have to be adjusted by chosing masses and couplings according to the relevant NC and 
CC values as listed in [I]]. For the CC case where Q\ should be identified with strange 

2 We could therefore, in principle, restrict the evolution of the charm FF to the nonsinglet sector. 


. The 

(18) 


9 








quarks the boundary condition in Eq. ( p~6|) is inadequate since m 2 a ~ A q CD is below the 
perturbative regime of QCD. We will have recourse to standard strange seas from the 
literature [|(| [27] instead. 


3 The Charm Fragmentation Function in SI DIS 


We will investigate the charm fragmentation function in CC and in NC SI DIS. In CC 
DIS the charm production mechanism is undebated since charm is dominantly produced 
by scattering on light strange quanta. Our reasoning will therefore be that D c is directly 
accessible in CC DIS at relatively low spacelike momentum transfer. An extracted D c can 
then be applied to NC DIS, where it might give insight into the details of the production 
dynamics ||. Also a test of the universality |28|] of the charm FF measured in CC DIS and 
e + e~ annihilation m would be an important issue directly related to the factorization 
theorems E3] of perturbative QCD (pQCD). 


All e c parameters discussed below refer to a Peterson type P3[ functional form given 


by 

D c (z) = IV jz [1 - z _1 - £ c /(l 
where N normalizes D c to f dzD c (z) = 1. 


-l 


( 20 ) 


3.1 CC DIS 

In CC DIS one does not expect to gain much insight into the charm production process 

since charm is dominantly produced in scattering events on strange quark^] in the nucleon, 

thereby permitting an experimental determination of the strange quark content of the 

nucleon BB0- On the other hand the well understood production mechanism makes 

a direct determination of the charm fragmentation function feasible by measuring the 

3 We assume a vanishing Cabibbo angle. Our results remain, however, unchanged if the CKM 
suppressed d —> c background is included. 
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energy spectrum of final state charm fragments. This is obvious in leading order accuracy 
whereQ 

da LO oc s(x)D c (z) (21) 

is directly proportional to D c . More precisely, it is not z but the closely related energy of 
the c —> /iv decay muon which can be observed in iron detectors i I, I- The smearing 
effects of the decay complicates the determination of the precise shape of D c but only 
weakly influences an extraction of (z) [|5j] which is valuable information if physically mo¬ 
tivated one-parametric ansatze [^, |3(J for D c are assumed. At NLO the production cross 
section is no longer of the simple factorized form of Eq. (^l|) [|l2j and double convolutions 
(symbol <g> below) of the form of Eqs. ®. o have to be considered. However, to a 
reasonable approximation 

d<JNLO = ([s®d,a s +g®dcrg\®D c )(xiZ,Q 2 ) 

= da LO K(x, z, Q 2 ) 

OC s(x) D c (z) K(x,z,Q 2 ) 

- s(x) V X>Q 2 [D c ](z) (22) 


holds also at NLO accuracy within experimental errors and for the limited kinematical 
range of present data on neutrinoproduction of charm. In Eq. (^2|) the approximate 
multiplicative factor T> absorbs the precise A'-factor K(x, z, Q 2 ) obtained from a full NLO 
QCD calculation |L2|. T> is not a simple universal fragmentation function but a nontrivial 
process-dependent functional which is, however, mainly sensitive on D c and shows little 
sensitivity on the exact parton distributions considered. The occurrence of x, Q 2 and z in 
Eq. (]22j) as indices and as a functional argument, respectively, reflects the fact that the 
dependence on x and Q 2 is much weaker than is on z. Eq. (^2]) tells us that s(y) fixes 
the normalization of da once K is known. On the other hand K (or V) can be computed 
|T2l| from D c with little sensitivity on s(y), such that s(y) and D c (z) decouple in the 
production dynamics and can be simultaneously extracted. This point can be clearly 
inferred from Fig. 1 where it is shown that the wide spread of CC charm production 
4 We will suppress some obvious scale-dependences in the following formulae and in their discussion. 
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predictions which were obtained in [|1^] using GRV94 [j26| and CTEQ4 |27]] strange seas 
can be brought into good agreement by a mere change of the normalization given by the 
ratio Sgrv(x)/ s cteq4(x)- The remaining difference is not within present experimental 
accuracy which can be inferred from the shaded band representing a parametrization 
T2| of CCFR data [|Gj]. High statistics neutrino data therefore offer an ideal scenario to 
measure D c complementary to an extraction from LEP data on e + e _ —> DX JlC 


32 


This has been Erst noted in [|28| where also a successful test of the universality of the 
charm FF has been performed. With new data [|], ^J. 7f[] at hand and with a sounder 
theoretical understanding |L2| of neutrinoproduction of charm it would be desirable to 
update the analysis in [^]. Nowadays one can in principle examine the possibility of a 
uniform renormalization group transformation from spacclikc momenta near above the 
charm mass (uN —> DX) to timelike momenta at the Z° peak (e + e _ —> DX). In 
E3 |3l| . |32|| charm fragmentation functions extracted from LEP data have been tested 
against px and p distributions measured in photoproduction at HERA. However, we 
believe that a comparison to z differential neutrinoproduction data is worthwhile beyond, 
since the latter measure directly the fragmentation spectrum whereas px and p shapes 
are rather indirectly influenced by the precise hardness of the FF via their normalization 
EMI- LInfortunately up to now no real production data are available but only strange 
sea extractions B § 0 resulting from an analysis of the former. We therefore strongly 
recommend that experiments publish real production data such that the above outlined 
program can be executed with rigour. Here we can only find a e c parameter which lies 
in the correct ball park and examine a few points which will become relevant for an 
extraction of D c once data will become available. 


An outstanding question is the possible effect of a finite strange mass on the full 
semi-inclusive charm production cross section |2T[ including 0(a]) quark scattering con¬ 
tributions. By comparing the thick and the thin solid curve in Fig. 2 (a) it is clear that 
the effect of a finite m s can be neglected even at low scales and for a maximally reasonable 


5 Data from NuTeV 
further analysis fMj. 


is to be expected in the near future; /A events observed at NOMAD await 
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value of m s = 500 MeV. For the larger scale of Fig. 2 (b) the effect of chosing a finite m s 
would be completely invisible. A further question which might influence the extraction 
of a universal FF from neutrinoproduction is the one of the scheme dependence in han¬ 
dling final state quasi-collincar logarithms In (Q 2 /m 2 c ). If these are subtracted from the 
coefficient functions as discussed in section |273|, the subtraction defines a running of the 
charm FF which becomes scale dependent.^ according to Eq. (jl7|). In Fig. 2 we examine 
such resummation effects for CCFR kinematics ||. We use the same Peterson FF with 
£ c = 0.06 |35|, |T^| once for a fixed order calculation [|TzJ (solid lines) and once as the non- 
perturbative part on the right hand side of the entire c D FF on the left hand side of 
Eq. ([17]) (dashed curves). We note that towards intermediate scales around Q 2 ~ 20GeV 2 
one begins to see the softening effects of the resummation which are enhanced as com¬ 
pared to FOPT. However, as one would expect at these scales, the resummation effects 
are moderate and could be compensated by only a slight shift of the e c parameter which is 
therefore, within experimental accuracy, insensitive to scheme transformations. We note 


that - as was already shown in |L2|] - according to Fig. 1 an e c of around 0.06 which 
we took from an older analysis in |J5| seems to reproduce the measured spectrum quite 
well. For (E u ) = 80GeV, ( Q 2 ) = 20GeV 2 a value e c ~ 0.06 gives an average {z) ~ 0.6 
consistent with ( z) = 0.68 ± 0.08 measured by CDHSW ||. In |]J[] a distinctly harder 
value of £ c ~ 0.02 was extracted from LEP data on e + e _ —> D*X. If the latter fit is 
evolved down to fixed target scales it is - even within the limited experimental accuracy 
- incompatible with the CCFR neutrino data represented in Fig. 1. From Fig. 2 it is clear 
that the difference cannot be attributed to a scheme dependence of the £ c parameter which 
is too small to explain the discrepancy. It would of course be interesting to know how 
much the above mentioned smearing effect of the c —> /iv decay might dilute the discrep¬ 
ancy. In any case, charm fragmentation at LEP has been measured by tagging on D*' s 
whereas neutrinoproduction experiments observe mainly IPs through their semileptonic 

6 Evolving the charm FF we adopt for consistency the evolution parameters m Cj b and of the 

CTEQ4(M) [^7j parton distribution functions and we use Qo = m c . 

' A similar e c has been obtained in |^l] in a related scheme. The e c value in [^2| has no connection 
to a massive calculation and cannot be compared to the values discussed here. 
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decay-channel (dimuon events). ARGUS |I6| and CLEO [[D]] data at y/s ~ lOGeV indeed 
show |38j a harder energy distribution of D *'s compared to D' s. It seems therefore to be 
possible within experimental accuracy to observe a nondegeneracy of the charm fragmen¬ 
tation functions into the lowest charmed pseudoscalar and vector mesons. We note that 
an e c value around 0.06 which is in agreement with neutrino data on D-production is also 
compatible with the D energy spectrum measured at ARGUS where the evolution may 
be performed either via FOPT using expressions in || or via a RG transformation along 
the lines of [|T^, [KJ. If forthcoming experimental analyses should confirm our findings 
the lower decade m c (~ lGeV) —*• ARGUS(lOGeV) may be added to the evolution path 
ARGUS(lOGeV) ->LEP(M Z ) paved in 0] for the charm FF. 


3.2 NC DIS 


The fragmentation function of charm quarks observed in NC DIS is of special interest 
since it allows for directly investigating |2], [3J the charm production mechanism which is a 
vividly discussed issue in pQCD phenomenology [24]. mm- Whereas for intrinsic heavy 
quarks one expects to observe a Peterson-like hard spectrum attributed to the dominance 
of the leading order quark scattering contribution, one expects a much softer spectrum for 
extrinsic GF since the gluon radiates the cc pair towards lower energies (^) during the hard 
production process before the nonperturbative hadronization takes place. Experimental 
analyses have been performed in ||, |J and the steep spectrum^] (best visible in Fig. 6 in 
||) of observed D-mesons together with the missing of a hard component at larger z seem 
to give clear evidence for the dominance of extrinsic GF over QS which was excluded at 
the 5% level 0. In a complete VFNS the QS ® component makes, however, just one part 
of the 0(a]) structure functions in Eqs. (|l8|) , (|l9|). Especially, there also exists a GF ^ 
component, albeit with the leading log part of it subtracted. Since the subtraction term in 
Eq. ( PI ) is proportional to D c and therefore only removes a hard component from GF^ 
one expects the rise towards lower z to survive the subtraction. Furthermore a perturba- 

8 The variable xd considered in 11 differs slightly from ^ in definition. The variables are, however, 
identical at the 2% level Idol. 
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tive evolution of the charm fragmentation function D c (z, Q 2 ) might soften somewhat the 
hard QS term. 


These expectations can be quantitatively confirmed in Fig. 3 where we show for HERA 
kinematics [§ the total (solid line) normalized 0{a\) production cross section for trans¬ 
verse virtual photons (F L = 0) on protons. We also show the individual components 
contributing to it: The processes Yd cc (dot-dashed) and Y c c 9 (dotted; inch 
virtual corrections) correspond to the GF W — SUB g and QS^ — SUB q — SUBn terms, 
respectively, subtracted at p = Q. They are not physically observable and only sensible 
if they are added to the QS ^ Born term (dashed) as in Eqs. (p~8|) , (p~9|) . We have per- 
turbatively resummed all logarithms of the charm mass via massless evolution equations 
starting at the charm mass [<5 0 = Q o = m c] and using the standard boundary conditions 
in Eqs. (|16|), ([lTj) for e c = 0.06. Finite charm mass effects on the subtracted QS^' con¬ 
tribution can be inferred by comparing the thick and the thin dotted curves, where the 
m c —> 0 limit has been taken for the latter. As has been theoretically anticipated in [|41[ 
the charm mass can be safely set to zero in QS d) and the involved convolutions in Eq. 
(|5|) may be replaced by the massless expressions in |8], which simplifies the numerics 
essentially and which we will therefore do for the /i = 2 m c curve in Fig. 4 below. As 
also stressed in j41| it is, however, essential to keep the charm mass finite in the GF d) 
contribution since m c tempers the strength of the z'~ l u-channel propagator singularity. 
Obviously the 0(a]) result of an ACOT based calculation deviates essentially from the 
naive Born term expectation and it seems by no means legitimate to treat GF^ as a 
higher order correction here. Contrary to the corresponding inclusive results in 0 and to 
the expectations in (T^j also the subtracted QS^ term is numerically significant in the 
semi-inclusive case considered here. In the light of the huge 0(a]) corrections it seems, 
however, undecidable as to whether include or omit QS ^ at 0(al) without knowing 
0(al) corrections within ACOT. 


In Fig. 4 the resummed result (VFNS: solid lines) can be compared to unsubtracted 
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0{a\)^ GF (1) (fixed order, dashed line). We show the total 0(a\) VFNS result for 
the choices ft — Q, 2 m c . The modest scale dependence arises exclusively through the 
QS d) term. For any of the other curves a variation of // is completely insignificant 
and we therefore only show n = Q. A full VFNS calculation (solid lines) seems hardly 
distinguishable from fixed order perturbation theory (dashed line) within experimental 
accuracy. The two approaches are even closer if one follows the suggestion in |14| and 
does not include (dot-dashed line) the subtracted QS^ term at the level of QS ® + 
GF^\ The data points in the figure correspond to HI measurements |2| of D Q (circles) 
and D* + (triangles; both including charge conjugation) spectra. The measurement is 
restricted to rjo < 1.5. Extrapolating to the full phase space gives rise to large acceptance 
corrections which are, however, quite uniform J40| over the kinematical range considered 
and therefore have a minor effect on the normalized spectrum. Since FOPT and ACOT 
based calculations are very close it seems improbable that an experimental discrimination 
between the two approaches will be possible. The born term in Fig. 3 is far from being 
the dominant contribution and an intrinsic c(x, Q 2 ) stemming from the resummation 
of perturbative logs can therefore not be excluded. The tendency of the data appears 
somewhat softer than any of the calculations and the tendency seems to be confirmed by 
preliminary data in a lower z bin |43|. The resummed calculation appears to be too hard 
at larger z around 0.6 if the QS^ component is included. As already mentioned, at the 
present stage of the calculations it can not be decided whether this hints at an intrinsic 
problem within VFNS calculations or wether this may be cured by 0(a^) corrections. 


4 Conclusions 


In this paper we have performed an ACOT \TQ based analysis of heavy quark fragmenta¬ 
tion in DIS including a calculation of semi-inclusive scattering on massive quarks at 0(a]) 
for general masses and couplings. As in the inclusive case [jlj effects from finite initial 
state quark masses can be neglected for practical applications to charm production in CC 


An 0(a^) NLO calculation h2[ within FOPT gives very similar results |40|. 
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and NC DIS. The involved convolutions in section [2.1| can therefore safely be replaced by 


their analogues in 


2^ , |l2|| . Neutrinoproduction is an ideal environment to extract the 
charm FF within DIS and a Peterson |23j type FF with e c ~ 0.06 seems to lie in the cor¬ 
rect ball park, where the sensitivity on the choice of scheme is small and finite m s effects 
are irrelevant. The e c value above is compatible with e + e - data if a nondegeneracy of 
charm quarks fragmenting into D's and D *’s is allowed for. For NC DIS it seems unlikely 
that a discrimination between fixed order and resummed calculations will be possible at 
HERA. Both approaches give similar results which show a spectrum that is somewhat 
harder than the tendency of the data [^. |3|. [43| . The resummed calculation is made worse 
if the 0(a]) quark scattering contribution is included at the perturbative level considered 
in this paper. 
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Appendix A: 

Matrix Elements for Real Gluon Emission 


The projections of the partonic Matrix Element onto the structure functions are most 
conveniently given in the Mandelstam variables A and which are defined below Eq. 
(0): 
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2(A 2 + m|S + _)' 

Sl 



(Al) 


where we conveniently use the shorthands A = A [m\, ml, —Q 2 ] and A' = A [ml, s, — Q 2 ]. 
In order to obtain the inclusive results f®(s i) in Appendix C of |I| the have to 

be integrated over 0 < y < 1, i.e. 

f?(si )=/ dy f?(s i,h) , (A2) 

Jo 

where y is dehned via the the partonic center of mass scattering angle 6* and related to 
t\ through 


V = ^ (1 + cosr) 

ii = > ^ 2 A' ( : y - y 0 ) , (A3) 

Si + 177-2 

with y 0 — [1 + (£ ++ + Si)/A']/2 being the would-be collinear pole of the t-channel prop¬ 
agator. 
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Figure Captions 


Fig. 1 The charm production cross section obtained in |T2| for two x, Q 2 points in the 
CCFR [6] kinematical regime. Up to a constant of normalization which was conve¬ 
niently chosen in Eq. (4) of [|l2j s e ff represents the triple differential cross section 
d 3 a/dxdydz where x and y are standard and z = Pd ■ Pn/q ■ Pn- Shown are the 
predictions using GRV94 [|6j] (solid) and CTEQ4 |27| (dashed) partons and a curve 
(dot-dashed) where the normalization of the CTEQ4 prediction is changed by mul¬ 
tiplying with the ratio of the strange seas sgrv(x)/ s cteq(x)- For all curves a 
scale-independent Peterson FF with e c = 0.06 has been used. 


Fig. 2 The same quantity as in Fig. 1. All predictions have been obtained for the CTEQ4 
parton distributions. The solid curves result from a scale-independent Peterson FF 
with £ c = 0.06. In (a) the thin solid curve has been obtained using the formulae 
of section [2] and chosing a finite strange quark mass of m s = 500MeV whereas the 
thicker curve corresponds to the asymptotic m s —> 0 limit (=MS). In (b) these two 
options would be completely indistinguishable and we only show the MS result. For 
the dot-dashed curves the final state quasi-collinear logarithm has been absorbed 
into a running of the charm FF. A Peterson FF with e c = 0.06 has been used as the 
nonperturbative part of the input as given in Eq. (|T7|). 

Fig. 3 The normalized charm production cross section da/dz (z = pn ■ Pn/q • Pn) for 
HERA kinematics (y/s = 314GeV). For a comparison with HI data p| the cross 
section has been integrated over lOGeV 2 < Q 2 < lOOGeV 2 and 0.01 < y < 0.7. 
Following the experimental analysis [2) only contributions from transverse photons 
{Fj j = 0) are considered. Shown is the total 0(a]) result (solid line) and the individ¬ 
ual contributions. Details to the calculation of the total result and the individual 
contributions are given in the text. The charm mass has been kept finite at the 
CTEQA value of m c = 1.6GeV everywhere except for the thin dotted curve where 
the m c —> 0 limit has been taken. 
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Fig. 4 A Comparison of the total 0{a\) result (solid lines) to HI data [Q on D° (circles) 
and D + * (triangles; both including charge conjugation) production. Shown are re¬ 
sults for two choices of the factorization/renormalization scale p. Also shown is the 
outcome of a fixed order 0{a\) GF calculation for comparison (dashed line). The 
dot-dashed line follows the suggestion in |I4|] and neglects quark initiated contribu¬ 
tions at 0(al), i.e. the difference between the solid (p = Q ) and the dot-dashed line 
is given by the (thick) dotted line in Fig. 3. For the dot-dashed curve as well as for 
the fixed order calculation (dashed) only p = Q is shown since the scale dependence 
is completely insignificant. 
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